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Spin-orbit coupled ultra-cold atoms provide an intriguing new avenue for the study of rich spin 
dynamics in superfluids. In this Letter, we investigate non-equilibrium spin dynamics in a spin- 
orbit coupled Bose-Einstein condensate (BEC) that are induced by sudden quantum quenches of the 
Hamiltonian. We observe a broad range of spin dynamics on different time scales, ranging from short 
time Zitterbewegung oscillations between two spin-orbit coupled bands, to long time spin relaxation 
and heating of the BEC induced by the interactions between atoms. We also demonstrate how 
quantum quenches can be exploited to populate a higher spin-orbit band, and observe a subsequent 
dipole motion. Our experimental results are corroborated by a theoretical and numerical analysis 
and showcase the great flexibility that ultracold atoms provide for investigating non-equilibrium 
spin dynamics. 
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The study of many-body dynamics far from equilib- 
rium has emerged as an important frontier in many 
branches of physics [1-7]. Ultra-cold atomic gases pro- 
vide a very promising setting for such studies because of 
the high tunability and precise characterization of sys- 
tem parameters [8, 9]. Recent experiments have opened 
a new avenue by demonstrating the realization of spin- 
orbit (SO) coupling in Bose-Einstein condensates (BECs) 
and degenerate Fermi gases (DFGs) [iqI-[i^. While 
SO coupling plays a prominent role in many important 
condensed-matter phenomena [l5l-[l7|. its realization in 
neutral atom superfluids is novel and provides a power- 
ful experimental platform for studying non-equilibrium 
spin dynamics due to a rich ground state phase diagram, 
intriguing spin dynamics and the presence of many-body 
interactions [18-31]. 

In this Letter, we investigate the non-equilibrium spin 
dynamics of a SO coupled BEC induced by quantum 
quenches of the Hamiltonian. The SO coupling is gen- 
erated by a Raman laser system, and quenches can be 
effected by sudden changes of Raman laser parameters. 
We observe short-time coherent spin dynamics as well 
as long time spin relaxation as a consequence of such 
quenches. Our main observations are the following: 

(I) On a short time scale (~ 500 /is), a quantum quench 
of the Raman detuning or a jump of the sign of the cou- 
phng strength induces spin oscillation between two SO 
bands, resembling the well-known Zitterbewegung (ZB) 
oscillation for relativistic Dirac electrons. The ZB os- 
cillation, first predicted by Schrodinger in 1930 |32j, de- 
scribes the fast oscillation or trembling motion of elec- 
trons arising from the interference between particle and 
hole components of Dirac spinors. Although fundamen- 
tally important, the ZB oscillation is difficult to observe 
in real particles and has only been measured recently 
using trapped ions as a quantum emulator of the Dirac 
equation [33]. In the SO coupled BEC, the two SO bands 
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FIG. 1: (Color Online), (a) Experimental configuration for 
the creation of spin-orbit coupling in the F — 1 manifold of a 
^^Rb BEC. (b) Typical band structure before (dashed-black) 
and after (solid-red) quenching the system by jumping the 
detuning S. The wavepackets symbolically show the wave- 
function directly after the quench. The short time dynamics 
are dominated by ZB oscillations, (c) Similar to (b) but for 
a jump of the relative phase between the two Raman beams. 
The band structure is unaltered by the phase jump. 



correspond to the particle and hole branches of the Dirac 
equation and the oscillation frequency is determined by 
the energy splitting between the two bands. 

(II) On the medium time scale the amplitude of the 
ZB oscillation decays away while atoms in the upper and 
lower bands move with different group velocities. The 
decay is accompanied by strong and sudden heating of 
the BEC and the many-body interactions between atoms 
reduce the decay of the oscillation amplitude. 

(III) On a long timescale (~ 20 ms) and in the con- 
tinued presence of evaporation from the dipole trap, the 
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heated cloud recondenses near the potential minimum of 
the lowest SO band, and the spin relaxes to the corre- 
sponding spin polarization. 

(IV) Furthermore, we show that quantum quenches 
can also be used to load a BEG into the upper SO band. 
We observe the subsequent dipole motion of the BEC 
as well as the accompanying change in spin composition. 
This measurement reflects the spinfullness of the upper 
band. 

Our experiments are conducted with BECs of ^''Rb 
of about 1 — 2 X 10^ atoms confined in a trapping po- 
tential with trapping frequency oJx^y^z = 27r x {20 — 
40, 174, 120} Hz, where the value of Ux depends on the 
intensity of the Raman beams (thus ft) as well as on 
a crossed dipole beam. Two crossed Raman lasers with 
wavelengths near A = 784 nm propagate along the e^^ie^ 
direction (relative angle = 90°), respectively. We apply 
a magnetic bias field of approximately 10 G in the ex di- 
rection (SO coupling direction) as shown in Fig. [TJi. The 
resulting quadratic Zeeman splitting for the F = 1 man- 
ifold Cz is 7.6£^r, which is sufficiently large such that the 
contribution of the hyperfine state |1,1) can mostly be 
neglected, yielding an effective spin- 1/2 system with the 
pseudo-spins defined as | t) = |1,0) and | |) = |1,— 1). 
Here Er = h^kl/2m = h x 27r x 1.866 kHz is the recoil 
energy and hkr = V27rh/X is the recoil momentum. 

In the pseudo-spin basis, the dynamics of the BEG can 
be theoretically described by an effective two-band ID 
Gross-Pitaevskii (G-P) equation with the corresponding 
Hamiltonian H = Hq -\- Vt -\- Hj. Here the single atom 
Hamiltonian is given by 
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after a local pseudo-spin rotation [10]. kx is the quasi- 
momentum of atoms, Q is the Raman coupling strength, 
and S is the detuning of the Raman transition from the 
level splitting, effectively acting like a Zeeman field. The 
harmonic trapping potential is Vt = ma;^x^/2. The 
many-body interactions between atoms are described by 
the nonlinear term 
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where the effective ID interaction parameters g^^ = 
and gn = flit = 9ii = with the spin- 

mayttz ^^'-^ ^'-^^ ^'-^'-^ mayttz ' ^ 

dependent 3D 5-wave scattering lengths cq and cq + C2 for 
Rb atoms (c2 = — 0.46ao and cq = 100.86ao). ao is Bohr 
radius, and z — \ — - — are the harmonic oscillator 
lengths. 

In our quench experiments we first prepare the system 
in the ground state with a finite detuning 5 and an initial 
quasi- momentum kx near kr- The system is quenched by 
either a jump of the Zeeman field from 5 to —5 (via a 
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FIG. 2: (color online) (a) Experimental realization of ZB oscil- 
lation for — 2.5Er and S jumped from 6A2Er to —6A2Er. 
Experimental data (open circles) is shown overlaid with re- 
sults of numerical simulations with no free fit parameters 
(black line) and analytic prediction based on effective two 
band Hamiltonian (dashed red line). Experimental (black 
dot) and numerical (solid line) quasi-momenta are shown in 
top part of plot, (b) Spin up population for numerical simu- 
lation of experimental parameters with (solid line) and with- 
out (dashed line) interactions, (c) Experimental image taken 
at t = 325/is showing the Stern-Gerlach separation and 2fikr 
photon momentum separation of the bare states, (d) Ex- 
perimental (open circles) and numerical results (filled circles) 
for ZB oscillation frequency vs. S after quench. Dashed line 
shows band splitting calculated from effective two band model 
for Q — 2.5Er. Experimental error bars are determined from 
fit uncertainties, shot to shot variations of kx, and calibration 
uncertainty of Q. 



frequency jump of the lasers), or by a phase jump of tt 
of the Raman field (which is equivalent to jumping Q to 
—Q in the Hamiltonian ([1])) (Figs. [Tl3,c). The jumps in 
S or the sign of Q are effected in less than 10 /is which 
is much shorter than the Rabi cycling period and any 
relevant system dynamics timescale. After the quench 
we allow the system to evolve for a given evolution time 
before starting the imaging procedure. The imaging pro- 
cedure consists of jumping off the Raman coupling and 
external confinement, allowing 11.5 ms time-of- flight in 
the presence of a Stern-Gerlach field that separates the 
bare states in the vertical direction, and imaging the bare 
spin states. The images are oriented such that the hori- 
zontal axis coincides with the direction of the momentum 
transfer of the Raman coupling. 

A concrete example of the short time dynamics follow- 
ing a jump of S from 6.42 Er to —6A2Er with Q = 2.5Er 
is shown in Fig. [2^. The dynamics are characterized by 
a rapid oscillation between the | |> and | t> spin states 
with a momentum transfer of 2hkr as seen in Fig. [2t. 
For the chosen parameters the oscillations occur at a fre- 
quency of 3.62Er and can be observed for many cycles. 
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The quasi-momentum remains relatively constant over 
this timescale. The role of interactions is revealed in the 
numerical simulations shown in Fig. [2]3, where it can be 
seen that the presence of repulsive interactions increases 
the duration of the oscillation The dependence of 

the spin oscillation frequency on the parameter S is plot- 
ted in Fig. [2]i along with the band excitation frequency 
and the results from numerical simulations of the non- 
linear G-P equation. Clearly the oscillation frequency is 
well described by the energy splitting between the lower 
and upper band. The occurrence of such oscillations is 
not unique to quenches of 5 and can also be observed 
under many other circumstances. For instance, we have 
observed similar spin oscillation for jumping 1] to —1] 
(see Fig. [5^). 

The spin oscillations can be understood theoretically 
as follows. Since after the quench the initial wavefunction 
is no longer an eigenstate of the new Hamiltonian, there 
is a fraction of BEC projected into the upper band. The 
components in the upper and lower bands beat against 
each other, leading to the observed oscillations in the 
bare state basis. It is straightforward to calculate the 
spin population N(j/ {N^-i-N^) = acr^ha- cos(A^t), where 
= Y^^(4^^^^^7^r^+~^P~+^^ is the energy difference of 
the two spin-orbit bands after the quench, the coefficients 
ttcr and ha- are determined by the system parameters [35] . 

The observed spin oscillation is a low-temperature ana- 
log to the well-known ZB oscillation studied in various 
systems 3^, 36|-41| because of the similarity between the 
effective two-band model in Eq. ([T]) and the Dirac-like 
equation. From the single particle Hamiltonian ([1]), the 
velocity operator of atoms in the Heisenberg picture is 
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FIG. 3: (Color Online). BEC loaded to Q = 1.8^r and 
5 — l.QEr^ followed by jump of 5 to —l.QEr. a) Experi- 
mentally observed quasi-momentum for the component of the 
BEC near kx — kr. b) Corresponding experimental (open 
circles) and numerical results (solid line) for spin composition 
of the BEC. Insets of a) and b) show images of the BEC at 
t = 8ms and t — 9ms (i.e. just before and after onset of 
heating). Vertical line indicates zero kinetic momentum, c-e) 
Numerical simulations show the real and momentum space 
composition for 2 ms , 6 ms, 8 ms, respectively. Black lines 
are the spin down component and red lines are the spin up 
component. 



Vx = K^x + krcrz)/m, (3) 

which does not commute with the single particle Hamil- 
tonian Hq and thus is not a constant. A straightfor- 
ward calculation [slj shows that the spin polarization 
{gz) actually oscillates with time, leading to a velocity 
oscillation (the quasi-momentum q = (kx) is roughly a 
constant on the short time scale) with a frequency de- 
termined by the energy splitting between two SO bands. 
Such a velocity oscillation yields the ZB trembling mo- 
tion of atoms. The amplitude of the position oscilla- 
tion of the BEC for our chosen parameters is small (be- 
low 1 /im) [35] and thus is not directly observable in 
our experiments. However, the oscillation of the velocity 
{vx) = ^{q^ kr {crz)) = {N^v^ + N^v^) /N is the veloc- 
ity of the center of mass of the BEC and is directly given 
by the spin polarization and which are measured in 
the time-of- flight images. Here the velocity oscillation 
amplitude ~ hkr/m ^ 4.14 mm/s. va^ = ^{q ^ kr) and 

= -^{q — kr) are the velocity of the two spin compo- 
nents in the time-of- flight. 

An ideal sinusoidal ZB oscillation is possible only when 
a single momentum is involved in the initial ground state. 



However, in a realistic system in a harmonic trap, the 
ground state has a small spread of the momentum around 
the minimum of the band, which we denote as g{k). 
The true dynamics of the ZB oscillation then become 
ilj{t) = Jj^ q(k)iljk{t)^ leading to a damping of the ZB os- 
cillation [40] on a short time scale (a few oscillation peri- 
ods). In our system, another source for damping comes 
from different group velocities of atoms in the two SO 
bands. The ZB oscillation is only present when the wave- 
functions in two SO bands have significant overlap in real 
and momentum space. When the wave packets in the 
two SO bands move with different velocities, they may 
start to separate in real and momentum space and the 
ZB oscillation starts to have strong decay ^] on a rel- 
atively long time scale. On the other hand, many-body 
effects expand the wavefunction in real space (increas- 
ing the wavepacket overlap) and narrow it in momentum 
space (reducing the momentum spread), leading to a re- 
duced damping effect on both short and long time scales 
as seen in Figs. [2)3 and [3)3. Therefore the inclusion of 
interactions into the system increases the duration of the 
ZB oscillations. Similar arguments apply to the ZB os- 
cillations generated by jumping the sign of the Raman 
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FIG. 4: Long time dynamics after a quantum quench. The 
BEG is loaded to ^ = 1.6Er and Q = lAEr, followed by 
a jump of 5 to —1.6Er. a) Quasi- momentum of the BEG 
after the quench (solid circles) , where the gray region indicates 
strong heating and rethermalization. For comparison, quasi- 
momentum of the BEG without quench is also shown (open 
triangles). Solid line indicates fit to dipole oscillation, b) Spin 
up population. Insets are experimental images taken during 
the evolution where the vertical line indicates zero kinetic 
momentum. 



strength Q (see Fig. [5^). 

After many milliseconds following the quench, the ZB 
oscillations decay while the BEG moves within the new 
band dispersion of the quenched Hamiltonian with the 
associated dipole oscillation in the lowest band. An in- 
teresting observation is the slow onset of dipole oscillation 
(~ 8 ms) and the subsequent rapid heating of the BEG, 
an example of which is shown in Figs. [3^ and [3]3. This 
data sequence was taken at 1^ = l.SEr and quenching 
S from l.QEr to —l.QEr. Note that the ZB oscillation, 
which typically damps out on the time scale of a few ms, 
is present in these cases as well but is not resolved in the 
experimental data shown in this figure due to the chosen 
1 ms time steps for this data sequence. The slow onset 
is also observed in the numerics and is related to the ge- 
ometry of the band structure. For a small value of the 
Raman coupling, the slope of the lower band after the 
quench of S is relatively flat. Thus the majority compo- 
nent of BEG in the lower band has a very small initial 
group velocity. We also verified in our numerics that for 
large Raman coupling, there is no such slow onset. 

The observed strong heating may also be explained 
from the geometry of band structure. The fractions of 
the BEG in the upper and lower band move with different 
velocities. When the difference in the group velocities 
are large, strong inelastic collisions occur and lead to the 
heating of the BEG. This physical picture is captured in 
our simulations (Figs. [St-e): at 6ms the different slopes 
of the two bands lead to a separation of wave packets in 
momentum space. Subsequently the wavepacket in the 
upper band turns around, and collisions with atoms in 
the lower band lead to the excitation of many momentum 
modes observed in the numerics 
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Once the BEG starts moving within the quenched band 
dispersion, it performs dipole oscillations 11, |43|. In 
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our experiments, the dipole oscillations occur as a con- 
sequence of a quantum quench, with initial wavepacket s 



FIG. 5: Loading of the upper band via two phase jumps 
separated by AO/us where Q = 3.5Er and 6 = l.GEr. 
a) Experimentally observed ZB oscillation (open circles) after 
a single phase jump, and corresponding numerical simulations 
(solid line) for experimental conditions, b) The time sequence, 
c-d) Experimentally observed quasi-momentum and spin po- 
larization, respectively. Insets are experimental images taken 
during the evolution where the vertical line indicates zero ki- 
netic momentum. 



in both lower and upper SO bands. An example of the 
dipole excitations generated by jump of 6 is shown in 
Figs.Sti andllD. Here the BEG is loaded with Q = lAEr 
and 5 = 1.6 Er and quenched to — 1.6 Er. Unlike dipole 
excitations well above AEr or in the upper band (Figs.[5t- 
d), here the BEG initially moves away from kx = pre- 
sumably due to momentum exchange with atoms in the 
spin up state (that are separated in momentum space by 
2hkr due to the Raman transfer), and becomes highly 
excited. In the presence of continued evaporation from 
the dipole trap, the cloud then recondenses near the new 
global minimum, with a small amplitude dipole excita- 
tion remaining as seen in the plot of the quasi-momentum 
in Fig. [4^. The frequency of the dipole oscillation around 
the new global minimum, at these parameters and for this 
small amplitude, is near the trapping frequency in the 
spin-orbit direction [11]. The spin relaxes to the almost 
fully spin-polarized state (Fig. HJd). 

In our discussion so far, we have concentrated on the 
dynamics following a single quantum quench. As we will 
show in the following, a sequence of jumps of the sign of 
Q can be used to load the BEG nearly entirely into the 
upper band and thus forms a pathway to studying upper 
band dynamics. The particular case that we choose to 
demonstrate this, is by two phase jumps of the Raman 
lasers: the first jump to start ZB oscillation, after which 
we wait for half a cycle until the vast majority of the pop- 
ulation has been transferred from | |) to | t) , followed by 
a second jump from —ft to Q. The corresponding time 
sequence is shown in Fig. [Sb. We experimentally inves- 
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tigate this procedure using Q = 3.5Er and 5 = l.dEr, 
inciting large ZB oscillation amplitudes by a phase jump 
of TT of the Raman fields, as seen in Fig. [5^. When the 
phase change is reversed by a second jump at t^=40 /is, 
we load approximately 80% of the BEC into the upper 
band near kx = kr. This is in agreement with an argu- 
ment based on the corresponding transformations in the 
Bloch sphere. As the population transfer is not unity 
for these parameters, there will be a residual ZB oscil- 
lation after the second jump. Allowing evolution time 
tw — tz^ the BEC undergoes dipole motion while it grad- 
ually melts. The quasi- momentum and the spin polariza- 
tion, defined as {N^—Ni)/{N^-\-Ni) for the component of 
the BEC in the upper band are plotted in Figs.[5fc and[5li. 
The change in spin polarization with quasi- momentum is 
consistent with the prediction from the a single particle 
model. For this calculation we fit the quasi-momentum 
measurement with a polynomial (shown as a line in Fig. 
[5fc). Using the single particle band structure we then 
calculate the corresponding spin polarization (shown as 
a line if Fig. [5]i). The good agreement with the experi- 
mental data indicates the spinfulness of the upper band. 

In summary, we observe rich non-equilibrium many- 
body spin dynamics in a spin-orbit coupled BEC after 
quantum quenches of Raman dressing parameters. The 
results presented in this work showcase the exceptional 
flexibility that cold atoms provide for the study of quan- 
tum dynamics in spin-orbit coupled superfluids. The rich 
physics accessible by rapid quenches of various system pa- 
rameters offers exciting outlooks for further studies, such 
as upper band dynamics, spin decoherence, etc. 
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Spin oscillation after the quench in Schrodinger picture 

For simplicity, in the following discussion we use the dimensionless single particle Hamiltonian, i.e., we measure 
energy in unit of Er, and kx in unit of kr. We assume that the Raman coupling strength is real. As an example, we 
consider quenching the detuning 5. Before the quench. 

The initial wave function is the ground state of this Hamiltonian, and kx is determined by the global minimum of 
the lower band dispersion relation. The ground state \ip{t = 0)) = \g) = (cos^ sin^)^, is a function of the system 
parameters Q and kx- After the quench, the Hamiltonian 

The eigenvalues are E± = /c^ + 1 ± ^ with = ^ {4kx — Sy + 11^, and the corresponding orthogonal eigenvectors 
are denoted as |+) = (cos a s'ma)^ and |— ) = (sino^ — coso^)^. Therefore we have (+|^) = (^|+) = cos(a — 0), 
{—\g) = (^'l— ) = sin(ce — 6). The time evolution of the wave function is 

\m) = e-'"'m = 0))=e-^"'\g) (6) 
= e-^*(l+)(+l + |-)(-|)l5> (7) 
= e-^^+*(+|5)|+)+e-^^-*(-|5)|-) (8) 

= cos(a -0) e-^+* + sin(a - 0) ( '^"^ ] e^^"* (9) 

^ ^ ysma J ^ ^ y— cos a J ^ ^ 



cos{a — 0) cos{a)e *^+^ + sin(ce — 6>) sin(a)e 
cos(ce — 0) sin(ce)e~*^+* — sin(ce — 0) cos(ce)e~*^-^ 



(10) 



Therefore the spin populations can be obtained as 



^ ^ = cos^(a-6>)cos^a + sin^(a-6>)sin^a+ ^sin(2a)sin(2a-26>)cos(A^t) (11) 

^ = cos^{a-0)sin^a-\-sin^{a-0)cos^a-^sm{2a)sm{2a-20) (12) 

and the spin polarization is 

P = -1 -± = cos{2a - 20) cos{2a) + sin(2Q^ - 20) sm{2a) cos{AEt). (13) 

Spin oscillation after the quench in Heisenberg picture 

The Hamiltonian after the quench can be written as 

H={kl + l)l2 + {2k:,-^)a, + ^a„ (14) 

where /2 is a 2 x 2 unit matrix. In Heisenberg picture, the wave function does not change while the spin operator 
evolves as 

^ = -iW^,H] (15) 

= -i{a,H-Ha,) (16) 
= -i{{a,,H}-2Ha,). (17) 

Using the fact that {gz^ 1} = ^cr^^, {az^ cr^} = 0, {<Jz^ cFz} = 2, we obtain 

'^=2i{{H-kl-l)cj,-{2k^-^-)}. 
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Integrating this differential equation, we obtain the formal solution for , 



IK 



+ 



2k, 



H-kl-l 



exp {2i{H -kl- l)t) 



(18) 



where (Jz{0) is the spin operator at t = 0. The last term is a fast osciUating function of time, leading to a fast 
oscillating velocity, according to Eq. (3) in the manusript. Such a velocity leads to the oscillation of the position, 
which is the Zitterbewegung oscillation. 



Using the fact that (+|a^(0)|+) = -(-|cr^(0)|-) = cos(2a), and (+|cr;,(0)|- 
calculate the expectation value of the spin polarization operator 

= cos2(a - e){+\cj,{t)\+) + sin2(a - e){-\a,{t)\-) 

+ isin(2a-2^){-|a^(t)|+) + isin(2a-2^){+|a,(t)|-) 



|c2(0)|+) = sin(2Q;), we can 

(19) 
(20) 
(21) 

(22) 



= cos^(a — 6) 



E+-kl-l 



2kx ^ 



E_-kl 



cos(2a) — 



■ cos(2ce) 



2 A^T" 



E+-kl-\) 
~ E_-kl-l 



^2i(E--kl-l)t 



+ i sin(2Q^ - 20) sm{2a)e'^'^^+-^--^^^ + ^ sin(2Q^ - 20) sm{2a)e'^'^^--^^--^^\ 



If we express the eigenvectors |+) and |— ) explicitly using the system parameters, we find that 

cos(2q^) 

and the spin polarization 



Oh — ^ OJc — ^ 

E+-ki-i ^ ~ E^-ki-r 



P = -J -± = cos{2a - 20) cos(2q^) + sin(2Q^ - 20) sm{2a) cos(AEt), 



which is exactly the same as the result in the Schrodinger picture. 
The ZB oscillation frequency for quenching S is (including unit) 



and the velocity oscillation term is 



huJzB = AE= y/{4.kxEr/kr-5)'^ ^Vt'^. 



. . , ,AE , 

— - sm(2a - 20) sm(2a) cos(— — t). 
m h 



So the maximum amplitude of the position oscillation is smaller than 

^ sin(2Q^ - 2(9) sin(2Q^) - —f^ — 

Trh hkr 
Er m 

- 1/kr, 



(23) 

(24) 
(25) 

(26) 

(27) 



(28) 



(29) 

(30) 
(31) 



which is much smaller than 1 fim and cannot be observed i n experiments directly. 



